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Abstract 

We present a complete analytical solution of a system of Potts spins on a ran- 
dom fc-regular graph in both the canonical and microcanonical ensembles, using the 
Large Deviation Cavity Method (LDCM). The solution is shown to be composed 
of three different branches, resulting in an non-concave entropy function. The ana- 
lytical solution is confirmed with numerical Metropolis and Creutz simulations and 
our results clearly demonstrate the presence of a region with negative specific heat 
and, consequently, ensemble inequivalence between the canonical and microcanoni- 
cal ensembles. 
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PACS numbers: 

05.20.-y Classical statistical mechanics. 

05. 70. -a Thermodynamics. 

89. 75. He Networks and genealogical trees. 



1 Introduction 

When a system phase-separates, it pays for tlie different domains with a surface 
energy, which is usually negligible with respect to the bulk energy. As a consequence, 
any non concave region in the entropy vs energy curve has to be replaced by a straight 
line. This is the result of the usual Maxwell construction. 
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However, the condition of negligible surface energy is violated in presence of long 
range interactions, as well as for systems with a small number of components. In both 
cases, the possibility of non concave entropies and ensemble inequivalence is well 
known, and has been demonstrated on numerous models, for instance [3,4,7,8]. The 
same condition of negligible surface energy is also violated on sparse random graphs: 
despite the fact that each site has only a small number of neighbors, there will be in 
general an extensive number of links between two (extensive) subsets of the system. 
The possibility of ensemble inequivalence in this type of models has been alluded to 
in some works related to the statistical physics of random graphs and combinatorial 
optimization [1]. However, these authors study the analog of the canonical ensemble, 
and replace the non concave part of the entropy by a straight line. This phenomenon 
remains thus to our knowledge unstudied, despite the widespread current interest 
in complex interaction structures, and networks in general. The purpose of this 
work is to present a simple, exactly solvable model on a random regular network, 
that displays a non concave entropy and ensemble inequivalence. This is a first step 
towards the study of more complicated networks, which may also include some local 
structure, like small world networks. 

The paper is organized as follows: in section 2, we present the model, and give its 
analytical solution; we then turn in section 3 to the comparison with microcanon- 
ical simulations using both Creutz [5] microcanonical dynamics and Metropolis [6] 
canonical simulations. The final section is devoted to conclusions and perspectives. 



2 Presentation of the model and analytical solution 

2.1 The model 



We study a ferromagnetic system of Potts spins with three possible states (a, b and 
c). The Hamiltonian is chosen to be: 

where (i, j) denotes all the bonds in the system, qi is the state of spin i, and 5q^q. 
is a Kronecker delta. In this form, the Hamiltonian simply counts the number of 
bonds between spins in different states. The ground state energy is 0. The spins are 
located on the nodes of a regular random graph where each node has connectivity fc, 
of order 1. A mean field like version of this model, with an all-to-all coupling, has 
been studied by Ispolatov and Cohen [7], and displays ensemble inequivalence. 
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Fig. 1. Schematic representation of the iteration (left), hnk addition (center) and site 
addition (right). Red nodes and sohd edges represent the original cavity spins and links, 
while the green colored nodes and dashed lines identify the additions. 



2.2 Analytical solution 



Random regular graphs possess very few loops of size order 1, and locally look like 
trees; this feature allows us to use standard statistical physics methods, originally 
developed for Bethe lattices. These calculations are usually done in the canonical 
ensemble only; in contrast, we are interested also in the microcanonical solution. 
We compute here the free energy and the entropy of the system, by following the 
formalism of the Large Deviation Cavity Method described by O. Rivoire in [2]. We 
consider however only large deviation functions with respect to spin disorder, and 
not with respect to disorder in the graph structure like in [2]. 

We call cavity sites sites which have only k — 1 neighbors, and one free link. Cavity 
site i sends a field hi along each link, which tells its state a, b or c. These field are 
distributed according to the probability distribution P {h): 



The first step is to obtain a self consistent equation for the probabilities p^, pb and 
through the analysis of the "iteration" process, represented on the left side of Fig 1. 
During an iteration step, a new site is connected to k — 1 cavity sites to become a 
new cavity site. Several possibilities must be accounted for, corresponding to all the 
possible configurations along the newly created edges. Let us note that for infinite 
temperature, or /3 = 0, each new spin has probability 1/3 to be in each of the three 
states a, b and c. This is the origin of the 1/3 factors in table 1 where we represent 
all the terms to be considered in the k = 3 case. 



P (h) = PaSh,a + PbSh,b + PcSh,c ■ 
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Table 1 

Analysis of the iteration process for k = 3: energy shifts and probabilities, ho is the field 
sent by the new cavity site. 



Using this table and following [2], we obtain: 
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from where we can easily calculate numerically Pa,b,c- ^oi larger k the generalization 
is straightforward, we have: 



Pa + {Pb + Pc) e 



k-1 



(3) 
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We compute the generalized free energy T {l3) through the formula: 



In 



+ - In 
2 



(4) 



where AEgHf. and AEnnk are the energy shifts due to a site and a link addition 
respectively. The ( . ) symbol denotes the expected value. Link and site additions 
are depicted on the center and right sides of Fig. 1, respectively. The analysis of the 
energy shifts in the k = 3 case is detailed in Tables. 2 and 3. 
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Table 2 

Configurations (/ii,/i2)) energy shifts IS.E and total probabilities Pi {AE) for the case of 
a link addition. The numeric factors stem from combinatoric arguments. 
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Table 3 

Possible configurations (/ii, /i2, /is), energy shifts AE and probabilities for the different 
states in which the new site can be. The overall factor of | corresponds to the a priori 
probability that the new site is in state a and the remaining numeric multipliers stem from 
combinatorics. 
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Fig. 2. Left: the three branches of the generahzed free energy ^ as a function of the 
inverse temperature /3, for k = 4. Right: the corresponding three branches for (3 (e) in the 
microcanonical ensemble. 



Plugging all the previous results in to Eq. 4, we obtain the expression of the gener- 
alized free energy of the system for the general k case: 
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where the three densities pa, Pb and pc are solutions of Eq. 3. Notice that this 
procedure does not necessarily yield a unique "free energy" T {(3)\ rather, there is 
one value of for each solution of the consistency equation (3). We must then 
follow all branches of the multi- valued function T (/3) to reconstruct the entropy 5* (e) 
through a generalized inverse Legendre transform (see for instance [9] for a use of 
this procedure in the context of signal processing): 



(6) 



where: 



e = 



'dp 



can easily be calculated numerically using finite differences. This is the final, implicit, 
solution for the entropy 5'(e). In fig. 2, we plot the different solution branches of 
and the inverse temperature /? (e). One clearly sees a negative specific heat 
region, signaled by the presence of multiple function values for the same energy. 
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3 Comparison with numerical simulations 

In this section we compare the analytical solution with the results obtained through 
numerical simulations. Microcanonical simulations were performed using Creutz [5] 
dynamics. During which, a fictitious "demon" is introduced, carrying an energy 
edemon- At each step, a spin fiip in the system is attempted, and the corresponding 
energy change SE is computed. If 6E < 0, the move is accepted; if SE > 0, the 
move is accepted only if edemon > SE. In both cases edemon is then updated so 
that the total energy E + edemon is kept constant; the energy of the system E is 
then constant up to a 0{1/N). For long run times, the demon's energy reaches an 
exponential distribution P{edemon = e) cx exp(— e/T^), from where one can compute 
the corresponding microcanonical temperature = 1/ (3rnu of our system: 

(7) 



Results of the Creutz dynamics are plotted on Fig. 3 and compared with the analyt- 
ical solution of the previous section. The agreement between the two is very good, 
with the f3 vs energy curve clearly showing a region of negative specific heat. 

Finally, we performed canonical Metropolis [6] simulations and calculated the average 
energy in the temperature range where our results predict ensemble inequivalence. 
As expected, the canonical caloric curve obeys Maxwell's construction and clearly 
"jumps over" the region where the specific heat is negative. 



4 Conclusion and perspectives 

We have presented a complete canonical and microcanonical solution of the 3-states 
Potts model on /c-regular random graphs, and shown that this toy model displays 
ensemble inequivalence. 

There is little doubt that this result should generically apply to models on different 
types of random graphs, such as Erdos-Renyi ones, among others. We also expect 
to observe ensemble inequivalence on small world networks, since in these systems, 
the presence of random long-range links should prevent the system from separating 
in two different phases. 

Beyond the inequivalence between microcanonical and canonical statistical ensem- 
ble, non concave large deviation functions should be expected for some properties 
on random graphs. Fig. 4 of [1] gives an example of this. The present work provides 
an example where the Large Deviation Cavity method allows to deal with such a 
situation, and to compute the non concave part of the large deviation function. 
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Fig. 3. Comparison for the caloric curve (3 (e) between the analytical solution (solid lines), 
the Creutz dynamics results (stars), and the Metropolis simulations (circles) for k = 4. 
The Creutz simulations were performed on networks with N = 40000 sites, for 10® "Creutz 
steps" , and the results were averaged over 20 network realizations. The Metropolis results 
were obtained using 50 different networks with = 10000 nodes, by performing 10^° 
Monte-Carlo steps. In both cases, the size of the error bars is comparable to the size of 
the symbols. 
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